Astronomy & Astrophysics manuscript no. landin 


© ESO 2008 


February 5, 2008 





Non-gray rotating stellar models and the evolutionary history of 

the Orion Nebular Cluster 

^ ; N.R. Landin^'^, P. Ventura^, F. D'Antona^, L.T.S. Mendes^'^ and L.P.R. Vaz^ 

O . 1 Depto. de Fi'sica, Universidade Federal de Minas Gerais, C.P.702, 30161-970 - Belo Horizonte, MG, Brazil; 
' e-mail: nlandin@fisica.ufmg.br, lpv@fisica.ufmg.br 
O i- Osservatorio Astronomico di Roma Via Frascati 33 00040 MontePorzio Catone - Italy 
O ' e-mail: ventura@mporzio . astro . it , dantona@mporzio . astro . it 

Depto. de Engenharia Eletronica, Universidade Federal de Minas Gerais, C.P.702, 30161-970 - Belo Horizonte, MG, Brazil; 
e-mail: luizt@cpdee.ufmg.br 

, Received . . . ; accepted . . . 
T— I ■ 

! ABSTRACT 

' Context. Rotational evolution in the pre-main sequence is described with new sets of pre-MS evolutionary tracks including 
rotation, non-gray boundary conditions (BCs) and either low (LCE) or high convection efficiency (HCE). 

Aims. Using observational data and our theoretical predictions, we aim at constraining (1) the differences obtained for the 
' rotational evolution of stars within the ONC by means of these different sets of new models; (2) the initial angular momentum 
■ of low mass stars, by means of their templates in the ONC. 
Methods. We discuss the reliability of current stellar models for the pre-MS. While the 2D radiation hydrodynamic simulations 
(-H _ predict HCE in pre-MS, semi-empirical calibrations either seem to require that convection is less efficient in pre-MS than in 
^ . the following MS phase (lithium depletion) or are still contradictory (binary masses). We derive stellar masses and ages for the 
I ■ ONC by using both LCE and HCE. 

O ^ Results. The resulting mass distribution for the bulk of the ONC population is in the range 0.2— 0.4M© for our new non-gray 
models and, as in previous analyse, in the range 0.1— O.3M0 for models having gray BCs. In agreement with Herbst et al. (2002) 
^ ; we find that a large percentage (~70%) of low-mass stars (M < O.SM© for LCE; M < O.SSM© for HCE) in the ONC appears to 
be fast rotators (P<4days). Three possibilities are open: 1) ^-^70% of the ONC low mass stars lose their disk at early evolutionary 
phases; 2)their "locking period" is shorter; 3) the period evolution is linked to a different morphology of the magnetic fields of 
the two groups of stars. We also estimate the range of initial angular momentum consistent with the observed periods. 
Conclusions. The comparisons made indicate that a second parameter is needed to describe convection in the pre-MS, possibly 
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' related to the structural effect of a dynamo magnetic field. 
Key words. Stars: evolution - Stars: interiors - Stars: rotation - Stars: Hertzsprung Russell and C-M diagrams 



1. Introduction 

We compute new sets of pre-MS models with the ATON 
code for stellar evolution, the version presented in Mendes 
et al. (^99 ), which includes stellar rotation according to 
the description by Endal & Sofia p976fl . updated for the 
present work to employ non-gray boundary conditions by 
Allard & Hauschildt l(TM7jl and AUard et al. §000). 

Our main goal is to improve our understanding of the 
appropriate physical constraints to be used for a general 
description of the evolution of stellar structure and its an- 
gular momentum with time. We are particularly interested 
on the choice of the stellar initial angular momentum and 
its variation with time, and, also, on the importance of 
convection efficiency during the pre-main sequence. To do 
this, we check our choices with respect to sets of relevant 
observations. 

Send offprint requests to: N.R. Landin 



Since the pioneering work by Henyey et al. H1955|l and 
Hayashi (jlftfiljl it is commonly accepted that pre-MS stars 
derive their luminosity by gravitational contraction, with 
the only exception of the short deuterium burning phase. 
Derivation of masses of young stellar associations has then 
generally been made by standard hydrostatic stellar mod- 
els including deuterium burning, an approach that we will 
use in this paper. This procedure rests on the assumption 
that neither the residual accretion after the protostellar 
phase nor the uncertainty in the zero point of ages affect 
the results in a strong way. 

Nevertheless, the theoretical description of moderately 
low and low mass objects is affected by the first princi- 
ple uncertainties in the description of some physical in- 
puts, in particular opacities, convection, equation of state 
(D Antona I1998|l and treatment of boundary conditions 
(BCs, Chabrier & Baraffe [l997f) . Since the low mass stars 
in the pre-MS are fully convective and over-adiabatic, any 
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change in the convective model substantially alters the 
location of the track in the theoretical HR plane. The 
use of a less efficient treatment of convection leads to 
larger temperature gradients, so that, for a given lumi- 
nosity, the structure readjusts on a more expanded con- 
figuration, with a consequent shift of the track to lower 
effective temperatures (Toff) (e.g. D'Antona & Mazzitelli 
HM^nWn D'Antona & Montalban lM)^ . The path fol- 
lowed by the theoretical pre-MS tracks on the HR diagram 
is also dependent on the boundary conditions used to fit 
the numerical integration of the structural equations of the 
interior with the atmosphere. The use of a non-gray atmo- 
spheric treatment shifts the tracks to cooler ToffS within 
an extended interval of masses and ages (Montalban et 
al. I2U(J4|I . The effect due to the non-grayness of the atmo- 
sphere is, in many cases, overwhelmed by the uncertainties 
related to the treatment of convection, that has a similar, 
but even stronger effect on the tracks, with only a few 
exceptions that are relevant for this work, as will be dis- 
cussed in W.'6\ 

Young stellar clusters provide a unique opportunity to 
test stellar pre-MS models. Many studies in the past have 
been focused on the Orion Nebula Cluster (ONC), be- 
cause it contains thousands of pre-MS objects. Hillenbrand 
H1997|l measured the V and I magnitudes (and colors) of 
^^900 stars and located them in the theoretical HR dia- 
gram by using bolometric corrections and taking into ac- 
count the interstellar extinction for the determination of 
colors and bolometric magnitudes, fn spite of the non- 
negligible uncertainties weighting on the derivation of 
these stellar parameters (Hillenbrand et al. I1997|l . such a 
work has been widely used, in connection with theoretical 
pre-MS tracks, to infer important information concerning 
the cluster itself. The mass and age distributions and the 
slope of the mass function can all be estimated by infer- 
ring the appropriate mass and age for each observed star. 
It is clear that the results obtained with this approach will 
depend to a certain extent on the set of tracks used to per- 
form the analysis, and on the physical inputs adopted to 
calculate the evolution. This is confirmed by the fact that, 
specifically for the ONC, studies that used different sets 
of evolutionary tracks reached significantly different con- 
clusions, particularly with respect to the mass and age 
distribution, and, more important, to the age spread and, 
consequently, to the evolutionary history of the star for- 
mation process within the cluster (Palla & Stabler T999). 

In the past few years, new detailed observational stud- 
ies of the ONC have been undertaken, focused on the 
rotational properties of the stars. Stassun et al. |1999) 
and Herbst et al. H2()()2(l measured the rotational periods 
of ~400 stars belonging to the ONC. All of them are in 
the Hillenbrand H1997|l sample, so they can be located on 
the HR diagram. More recently, Stassun et al. 1)200411 and 
Flaccomio et al. I)2()()3all2008b|l reanalyzed all the archival 
Chandra observations of the ONC studying in great detail 
the X-ray properties of the observed objects, in an at- 
tempt to elucidate the origin of X-ray emission in pre-MS 
stars. All this available information renders the ONC an 



excellent laboratory to test stellar evolution theories of the 
pre-MS phase. 

None of the previous analyses of the ONC rotational 
database has been done using non-gray models, and the 
effect of using different convection efficiencies has not 
been extensively tested. We therefore decided to use this 
database to test and calibrate our new sets of non-gray 
tracks for rotating stellar models, with masses in the range 
0. 085 <M/M0< 3. 8. As a byproduct, we will have insight 
on how much the rotational properties of this popula- 
tion depend on the choice of the evolutionary tracks. At 
present, there is no definitive observational constraint that 
can be used to choose among the tracks obtained by using 
different convection inputs. In this work we derive masses 
and ages with different sets, in order to appreciate how 
much the results we are presently interested in, namely 
those on the angular momentum evolution, depend on the 
choice. In particular, we will test the role of (a) boundary 
conditions (gray and non-gray), (b) different convection 
efficiencies in the framework of the MLT, and (c) rotating 
and non-rotating models. 

As indicated by the non-negligible differences among 
the set of tracks adopted, we show that the detailed period 
distribution as a function of mass and age is dependent on 
the physical inputs. However, the qualitative information 
on the rotational distribution of stars of different mass in 
the ONC remains similar, and we confirm that the distri- 
bution of periods is bimodal only for masses larger than a 
"transition" mass depending on the convection model. 

In 321 we summarize the uncertainties in the theoretical 
models, and also try to clarify some issues left ambiguous 
in the recent literature on pre-MS evolution. In |3|we pro- 
vide a brief description of the latest rotational, non-gray 
version of the ATON stellar evolution code, and describe 
the choice of the main physical and chemical inputs used 
for the present work in ^ The role played by the convec- 
tive modeling on the location of pre-MS tracks on the HR 
diagram, the role played by non-grayness, and the lithium 
pre-MS depletion are discussed in iJHl In S|ni we present 
the photometric and rotational data related to the ONC 
stellar population, and the derivation of masses and ages 
from our new sets of tracks. The results on the rotation 
period distribution and its possible interpretation(s) are 
discussed in Conclusions are given in iJHl 

2. Uncertainties in the standard hydrostatic 
stellar models 

We briefly rewiew the reliability of the procedure of deriv- 
ing ages and masses of young stellar population based on 
hydrostatic evolutionary tracks. We first discuss two pre- 
liminary important physical points: (1) Is it appropriate to 
neglect residual accretion from the disk? (2) As the models 
disregard the hydrodynamic formation process, how do we 
take into account the zero point of ages? We summarize 
the well known uncertainties in the location of evolution- 
ary tracks in the HR diagram and show why current data 
do not allow us to choose among models. 
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2.1. The role of residual accretion 

The pre-MS phase begins after the main accretion phase, 
during which the protostellar core is formed and the star 
is embedded into the dust of the forming cloud. When ob- 
jects of low mass become luminous in the visible or near 
infrared bands, the main accretion phase can be consid- 
ered finished, although the accretion disk is still present in 
many cases. Accretion then will not change the final stel- 
lar mass much, while the presence of the accretion disk 
may be very important in determining the rotational evo- 
lution of the star. In the pre-MS, accretion still may occur 
at rates in the range '^3x10^^ — 4xlO^^MQ/yr (Basri 
& Bertout, [TU5^ for classical T Tauri (TT) star, and 
mostly below ~lO~"'^°M0/yr for very low mass stars be- 
low -0.2Mo (Mohanty et al. EOnS)- Walter l(TU57jl and 
Strom et al. I)1H8H(I noticed the coexistence of classical 
and weak TT stars in the same region of the HR diagram, 
corresponding to an age range (derived from standard hy- 
drostatic models without accretion) from ~3xl0^yr to 
~3xl0^yr, concluding that the disk evolution is decou- 
pled from the evolution of the central star, and that the 
disk may disappear either early or at late stages. This will 
affect the stellar rotation in different ways, as shown by 
subsequent studies (e.g. Bouvier et al. 11993(1 . interpreted 
in terms of magnetic disk locking (Bouvier et al. I1997|l . 

Can the residual mass accretion alter the stellar evo- 
lution? Deuterium burning from the accreted matter pro- 
vides a luminosity L--15L0X (Af/(lO^^M0/yr)) (Stabler 
I1988|l . which could inffuence the evolution only for very 
high rates of mass accretion, or at very low luminos- 
ity. Accretion itself may alter the evolution if the ac- 
cretion timescale M/Mc^t is of the order of the ther- 
mal (Kelvin-Helmholtz) timescale tKH- At high lumi- 
nosity, where txH^lO^yr for a typical TT, this con- 
straint is not respected only for very high accretion rates. 
At about the solar luminosity tKH"^2x lO^yr, for which 
M:rit^5xlO~^MQ/yr, which is at the upper boundary of 
the values observed in TT. In this case, or in similar cases, 
the evolution can be altered, but this ill cannot be sus- 
tained for a very long time, otherwise we would not find 
low masses surviving. Thus we conclude that evolution in 
which accretion plays a dominant role, such as described 
e.g. by Hartmann et al. 111997') and Tout et al. 1)19991 1 may 
be relevant only for a short lifetime, in a small fraction of 
the pre-MS stars, and that the bulk properties of young 
stellar associations can be derived by adopting traditional 
pre-MS hydrostatic models. 

2.2. The zero point of stellar ages 

One other item that must be taken into account is the def- 
inition of the zero point of ages, which, for the hydrostatic 
pre-MS evolution, is connected to the location in luminos- 
ity of the "starting model" , i.e. the internal thermody- 
namic conditions within the star when the mass accretion 
process ends. The uncertainty in the starting models is 
reflected in the age-luminosity relationship. A first uncer- 



tainty is due to the fact that we do not know whether the 
deuterium burning phase occurs during the visible pre- 
MS or during the main accretion phase: this depends on 
the average protostellar accretion rate (see D'Antona & 
Mazzitelli 119971 for a discussion). For the stars that do 
not completely burn deuterium in the protostellar phase 
(M ^ O.5M0), the typical uncertainty in the lifetime, how- 
ever, must be of the order of txH at the end of the main 
accretion phase (e.g. Tout et al. I1999|l . In the model by 
Palla & Stabler (il990, 199^, in which the stellar core of 
low mass stars emerges at visible wavelengths close to the 
deuterium burning region in the HR diagram (birthline) 
due to the "thermostat" action of D-burning, the typical 
txH is the thermal timescale at the D-burning luminosity 
and radius, ~10^yr for masses of a few tenths of M0, and 
up to ~10^yr for O.IM©. For the lowest masses, however, 
the "birthline" concept may not be valid, as the main ac- 
cretion phase may end before deuterium is ignited, and the 
stars probably start their hydrostatic evolution at higher 
luminosities, where the thermal timescale, and thus the 
uncertainty in the zero point of ages, would be smaller. 

On the other hand, we must be careful not to confuse 
the "numerical" uncertainties with the "physical" uncer- 
tainties. Hydrostatic contraction will let the star ignite D- 
burning as soon as the central physical conditions allow 
it, and it seems very unlikely that we can take as a zero 
point of the pre-MS a model with central temperatures in 
the middle of the deuterium burning phase, disregarding 
the action of the thermostat^. 

2.3. The role of non-gray atmospheres 

In the gray atmosphere approximation the increase of the 
temperature T from the surface moving inwards is de- 
scribed via a relationship between T and the optical depth 
(r), and the pressure is calculated by integrating the hy- 
drostatic equilibrium equation. The match between the 
interior and the external layers is made at t=2/3. Any 
link between pressure and temperature in the atmosphere, 
fixed by the either a convective or radiative gradient, is ig- 
nored. In the non-gray treatment a self-consistent integra- 
tion is performed down to an optical depth at which the 
diffusion approximation is valid (Morel et al. 11994,) . and 

^ Some models of the widely adopted grids by Baraffe et al. 
P 998(1 begin their evolution at a luminosity in the middle of 
D-burning, which lasts for ~2xl0^yr at O.IM0. This is due 
to the fact that the lowest available gravity in the NextGen 
(Allard & Hauschildt 119971 AH) atmospheric grid, which these 
models employ, was logg =3.5. The more recent models by 
Baraffe et al. P002,) start at lower gravity, and reach the D- 
burning luminosity by contracting in thermal equilibrium. Of 
course, in these two sets of models the temporal evolution of 
luminosity is very different for the first ~2xl0®yr, but this 
does not automatically mean that this figure is the physical 
uncertainty in the ages, as the "numerical" choice of starting 
from a model with central temperature in the middle of the D- 
burning stage does not necessarily correspond to a physically 
plausible behavior (Montalban & D' Antona l20(J6l . 
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including the treatment of atmospheric convection, which 
cannot be neglected at low ToffS. The use of frequency- 
dependent opacities may also modify the onset of con- 
vection within the atmosphere, and both the Tcs and the 
colors of the tracks can be strongly affected. The necessity 
of adopting outer boundary conditions based on realistic 
non-gray atmosphere models for the pre-MS and low mass 
MS was pointed out by Chabrier &; Baraffe H1997I and ref- 
erences therein) , who have shown that the use of radiative 
T(r) relations or gray atmosphere models is invalid when 
molecules form near the photosphere, at Tgg below 4000K. 
Outer boundary conditions based on the gray assumptions 
yield hotter models for a given mass. According to Baraffe 
et al. (^998) the use of an inappropriate outer boundary 
condition, such as the Eddington approximation, yields an 
overestimation of Tes for a given mass up to 300 K. 

2.4. The role of convection 

There is broad consensus in the literature that the treat- 
ment of superadiabatic convection in the pre-MS affects 
the tracks location to a great extent. In a series of works fo- 
cusing on understanding the role played by different physi- 
cal inputs on the pre-MS evolution, D'Antona & Mazzitelli 
ifTMTji . D'Antona & Montalban f)2003jl . and Montalban et 
al. H2()()4|l outlined the major impact of convection model- 
ing on the location of the stellar tracks in the HR diagram. 
Convection was found to be by far the most relevant ingre- 
dient influencing the determination of the mass and age 
of observed stars. 

When modeling convection, it is essential to specify the 
mixing scale A, i.e. the typical distance that convective 
eddies travel before dissolving and delivering their excess 
gravo-thermal heat to the environment. The role played 
by A is relevant for the determination of the temperature 
gradient, as the conservation of flux implies that a larger A 
must be compensated by a lower degree of overadiabatic- 
ity. The convective flux behaves as Fc^A'^ in zones where 
most of the energy is carried by convection. Conversely, 
within low efficiency convective regions, Fc^A^. The im- 
pact of the choice of A is therefore more evident where 
convection is not efficient. The pre-MS tracks result to 
be particularly sensitive to convective modeling, because 
the surface convection extends to most of the (if not the 
whole) star, and the low densities involved (particularly 
in the early phases of gravitational contraction) make the 
convective process highly inefficient. 

Presently, the main ways of computing convection in 
stellar envelopes, for wide grids of stellar models, are: 

1. The traditional mixing length theory — MLT (Bohm- 
Vitense 119581 and subsequent variations of this same 
model) — assumes that both the dimension of the 
convective eddies and the mixing length are propor- 
tional to the local value of the pressure scale height, 
i.e. A = aHp, where a is a free parameter that is usu- 
ally calibrated in order to reproduce the solar radius. 



2. In the Full Spectrum of Turbulence model — FST 
(Canuto et al. I1996|l — the whole spectrum of ed- 
dies' dimensions is considered, and the mixing length is 
taken as the distance of the nearest convective border. 

3. MLT, in which the a value for each gravity and T^s is 
calibrated upon 2D or 3D hydrodynamical simulations. 

Ludwig et al. (|1999;i . using their 2D radiation hydrody- 
namic models, have provided a calibration of the param- 
eter a in a wide region of Toff's and gravities^ to be used 
in the computation of gray stellar models. These 2D mod- 
els indicate that convection in the pre-MS is on average 
more 'efficient' than in the MS, corresponding to a larger 
a. The idea of calibrating the average a using numerical 
simulations has been extended by now to a few 3D com- 
putations: Ludwig et al. (AHS. 000211 . for an M dwarf at 
Teff = 2800 K and logg = 5, find a ~2.1; and Trampedach 
et al. H1999|l , for the range of main sequence gravities and 
log Teff =3. 68-3. 83, find a ~1.6-1.8 in the whole range. 
Asplund et al. H2000fl have compared 2D and 3D atmo- 
sphere models for Sun, and found that the 2D solar model 
has marginally larger gradients than the 3D one. Although 
an extrapolation to regions not explicitly computed is 
not allowed, also these few 3D models indicate efficient 
convection in the overadiabatic envelope. Montalban & 
D'Antona 1)2006(1 have computed gray models by using 
this calibration: the tracks they obtain are very similar 
to the FST tracks'^. Unfortunately, a very efficient convec- 
tion in pre-MS is not consistent with the lithium depletion 
patterns of young open clusters (D'Antona & Montalban 
I2003f) . Any attempt to calibrate convection efficiency in 
pre-MS by means of comparisons between binary masses 
dynamically determined and those assigned from different 
sets of evolutionary tracks (HCE and LCE) , seems to be 
ambiguous (Landin et al. 12006(1 . 

2.5. The role of convection coupled with non-gray 
atmospheres 

We note a point that is generally overlooked, but that is 
crucial in order to understand the relevant parameters in 
track building: Montalban et al. 1(2004(1 have shown that 

^ Q is mapped in the domain T^s = 4300 — 7100 K, logg = 
2.54 - 4.74 

^ The FST convection model corresponds to very efficient 
convection, as shown by the quasi-coincidence of the result- 
ing tracks with the tracks employing the MLT a calibrated 
on the 2D hydrodynamic models. We find confusing and mis- 
leading the statement by Baraffe et al. (I2UU2II . who point out 
that the FST model is very inefficient in the upper solar lay- 
ers, where it provides results that are not consistent with the 
hydrodynamic simulations of convection for the solar model 
(nevertheless, remember that the FST model provides a bet- 
ter fit than MLT for the spectrum of solar oscillations, see e.g. 
Canuto & Christensen-Dalsgaard, 19981. What matters in the 
description of the pre-MS is the average efficiency of convec- 
tion in the whole superadiabatic envelope, and this is very large 
for the FST, roughly corresponding to an a value in the MLT 
description somewhat larger than 2. 
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the problem of convection is intertwined in a subtle way 
with the problem of non-gray boundary conditions. For 
example, in the MLT framework, convection in the pre- 
MS should be described not only by the ratio a — I /Hp 
in the interior of the star (ain), below the non-gray atmo- 
sphere, but also by the value that this parameter has in 
the atmosphere itself (ttatm), and by the matching point 
between atmosphere and interior (Tph). In particular, the 
two widely used sets of model atmosphere by Baraffe et 
al. (1998| are referred to as the set built with a=1.9 and 
the set having q;=1.0. Thus they are supposed to pro- 
vide a clue to how the tracks vary by changing from a 
moderately high convective efficiency (HCE) represented 
by the set a=1.9, to a low convection efficiency (LCE) 
represented by the set a=l. However, in Baraffe et al. 
H1998|l . the parameter a refers only to the value of ain, 
and is misleading for two reasons: (1) the set aiii=1.9 
stops at masses M>O.6M0, and only the set Q;in=1.0 is 
available for masses M<O.6M0; (2) for both sets, the at- 
mospheric model grid adopted is the same, and computed 
with Qfatm = 1-0. Montalban et al. H2004(l have shown 
that the fact that most of the superadiabatic part of the 
envelope is computed with a very inefficient convection 
(aatm=l-0) shifts the Toff by ~ — 150K for the solar pre- 
MS. Thus one may be led to attribute the smaller TcffS of 
the Baraffe et al. (|1998|) tracks to the use of the non-gray 
atmospheres, whereas they are due in part to the fact that 
these non-gray atmospheres are computed with LCE. 

3. Rotating stellar models 

The pre-MS tracks were calculated by means of the ATON 
stellar evolution code. A full and detailed description of 
the numerical structure can be found in Mazzitelli ( 1989} 
and Mazzitelli et al. ( T99F) . Here we briefly recall the main 
micro- and macro-physics input. 

3.1. Rotation 

As described in Mendes et al. H1999|l rotation was im- 
plemented in the ATON code according to the ap- 
proach followed by Endal & Sofia H197fi|l . which uses the 
Kippenhahn & Thomas H197U|I method improved with a 
potential function that includes a term related to the dis- 
tortion of the figure of the star. In this approach the spher- 
ical surfaces used in standard stellar models are replaced 
by the equipotential surfaces. 

The current version of the code allows one to choose 
among three rotational schemes: 

— Rigid body rotation throughout the whole star 

— Local conservation of angular momentum in the whole 
star 

— Local conservation of angular momentum in radiative 
regions plus rigid body rotation in convective zones 

It is possible to consider only the hydrostatic effects of ro- 
tation and/or to include also the internal angular momen- 
tum redistribution and surface angular momentum loss. 



3.2. Atmospheric treatments 

In addition to the gray atmospheric boundary conditions 
in the version by Mendes et al. (1999J, the ATON2.4 
code can now adopt non-gray atmospheric integration. In 
the first case the internal structure is matched at t=2/3 
with the values of P and T found via a Krishna-Swamy 
((1966) T{t) relation. In the non-gray case we can choose 
the optical depth at which the matching is done: typical 
values are in the range l<r<100. For the atmospheric 
structure we can choose among several sets of models: 
(1) either the MLT or the EST models by Heiter et al. 
(2002); (2) the AH-hAHS grid, based on the NextGen 
models by Allard & Hauschildt (19!2|), complemented by 
the low gravity models by Allard et al. I)2000|l . To follow 
the evolutions starting from early, low-gravity stages we 
merged the NextGen tables (only available for log(?>3.5) 
with the more recent low-gravity models (2<log (7<3.5). 
Eventually, we obtained rectangular tables in the range 
2<log5<6 and 2000 K<roff<6800 K. 

3.3. Input micro-physics 

The radiative opacities are taken from Iglesias & Rogers 
(IR93, '19931, extended by the Alexander & Ferguson 
(AF94, ), 1994,1 tables in the low-temperature regime. The 
OPAL equation of state (Rogers et al. 119961 R96) is used 
in the range 3.7<logr<8.7, while in the low-T high den- 
sity regime we use the Mihalas et al. (M88, rr9B8jl EOS. 
The nuclear network includes 14 elements and 22 reac- 
tions; the relevant cross-sections are taken from Caughlan 
& Fowler 1,1988^) . 

3.4. Convection 

The borders of the instability regions are found via the 
Schwarzschild criterion. The convective fiuxes can be com- 
puted either by the EST model (Canuto et al. ^996), or 
by the MLT fBohm- Vitense [T^55|l : in the latter case the 
free parameter a = K/ Hp determining the mixing length 
can be arbitrarily selected. 



4. Input of present models 

We computed pre-MS stellar evolutionary tracks in the 
mass range 0.085<M/Afo<3.8. We adopted the solar 
chemistry with Z=0.0175 and F=0.27, while the starting 
deuterium abundance in mass fraction, following Linsky 
(1159 8), was set to X{D)^2xlQ-^. The main physical pa- 
rameters of the model sets are summarized in Tabled The 
evolutions start from a fully convective configuration with 
central temperatures in the range 5.3<logTc<5.8, follows 
deuterium and lithium burning and ends at the main se- 
quence configuration. 
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Fig. 1. Paths followed by the theoretical tracks (left) and isochrones (right) of the rotating models calculated with 
three different values of the free parameter a determining the mixing length. For clarity, we report only the following 
masses calculated with our three sets of tracks (al.0,a2.0,a2.2): 0.09, 0.1, 0.2, 0.3, 0.35, 0.5, 0.7, 1.0, 1.4, 2.0, 2.5, 3.0, 
3.3 and 3.8 Mq from bottom to top. The observational data from Hillenbrand H1997|l are represented by • symbols. 



4.1. Boundary conditions 

We compute gray models for comparison with previous 
results. For the non-gray models we use the AH-I-AHS 
grid. Following a suggestion by Heiter et al. (|2()02|l . we 
match the atmospheric grid with the interior integra- 
tion at rph=10. This choice should minimize the consis- 
tency problems related to the different EOS and opacities 
adopted in the interior and in the atmosphere, and to the 
absence of turbulence pressure in the atmospheric model- 
ing. 

4.2. Choices about convection 

As the FST atmospheric tables by Heiter et al. i2002p 
are available only for temperatures above rcff=4000 K, 
they cannot be used to compute models below ^ O.TM© 
(Montalban et al. I2004|l . In particular, they cannot be 

Table 1. Main physical parameters of the present models. 



Parameter 


Input 


Mass 


0.085-3.8 Mq 


Boundary Conditions 


AH4-AHS 


Matching point 


r=10 


Convection model 


MLT 


MLT parameter 


a=1.0, 2.0 and 2.2 


Rotation 


Rigid body 


Initial angular momentum 


Kawaler llTWt 


Opacities 


IR93 and AF94 


Equation of state 


R96 and M88 


Chemistry 


Y=0.27, Z=0.0175 


Initial X(D) 


2x10"^ 



used to analyze the ONC stellar population, mainly con- 
centrated at Toff <4000 K. To achieve consistency between 
the internal and the atmospheric convective treatment, we 
decided to calculate only MLT models. 

Apart from the quoted cases in which a can be cali- 
brated on hydrodynamic models, the necessity to simplify 
the numerical treatment of convection leads stellar mod- 
elists to use a single a for the convection zone and for 
all the evolutionary phases. This choice is equivalent to 
adopting an average efficiency of the convective transport 
on the whole extension of the convective region, and on 
all the evolutionary phases. There is no good reason to 
assume that the a that, e.g., fits the solar radius should 
be used for other masses and for different evolutionary 
phases. Further, the efficiency of convection might change 
considerably within a convective zone, thus requiring the 
use of a variable a. Thus a preliminary investigation of 
the effects of changing this parameter is mandatory. 

We computed sets of models with three different values 
of the parameter: the models with ain=2.0 (ck2.0 set) allow 
a fit of the solar radius for non-rotating models^; the mod- 
els with Q:in=2.2 (a2.2 set) are chosen to provide a "very 
efficient convection" set. Both these sets are termed HCE 
(high efficiency convection) sets. We further provide tracks 
with Q!in=1.0 (al.O set). This latter choice, according to 
D'Antona & Montalban (!2003', confirmed in ii5.4(l . leads 
to a better agreement with the lithium vs. Tag relation 
observed in young open clusters stars. Remember, how- 
ever, that the adopted BCs come from model atmospheres 
computed, down to Tph=10, with aatm=l-0. (see the dis- 

^ Also Baraffe et al. (fTM^ and Montalban et al. ipiTOjl find 
a similar ain (=1.9) to reproduce the solar radius. 
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cussion in §2.5 and in Montalban et al. I2(JU4|I . We show 
in Fig. ^ the comparison between tracks and isochrones 
of these three sets, showing the well known fact that both 
masses and ages are affected by the choice of the convec- 
tion model. 

4.3. Initial angular momentum 

Rotation was modeled according to the rigid body law. 
This choice is motivated by the fact that most of the 
low mass stars are still fully convective in the evolution- 
ary stages of interest for this study. As a first attempt, 
the initial angular momentum Jin for all models was esti- 
mated according to the prescriptions by Kawaler 1)1987(1 . 
In that work, a relationship between angular momentum 
and stellar mass for stars earlier than FO was derived us- 
ing main-sequence (MS) stellar models and an estimate 
of the mean initial angular momentum-mass relation was 
made for stars of later spectral type. As masses larger 
than about I.BMq do not lose much angular momentum 
during their early evolutionary phases, it can be assumed 
that these stars reach the MS with the same angular mo- 
mentum that they had at the beginning of their evolution. 
Kawaler H1987|l is able to reproduce the observational re- 
lation at M>1.5M0 using his own models for radiative 
stars, he then extends the models to lower mass stars. For 
these, however, the observations do not provide a direct 
comparison with the initial angular momentum, as their 
rotation has been slowed down during the main sequence 
lifetime (but not significantly during the pre-MS, if we 
adopt the hypothesis of pre-MS disk locking in the same 
way as Bouvier et al. ll997(l . For the range 0.6- 1.25M0, the 
initial angular momentum-mass relation can be easily ob- 
tained from the respective mass-radius and mass-moment 
of inertia relations from Kawaler H1987|l : 

Jkaw = 1.566 X 10^° f — j cgs. (1) 

Once started with this initial angular momentum, we 
keep it constant in our models during the pre-MS accord- 
ing to the above mentioned hypothesis, whose validity will 
be checked later. This expression will then be extended 
to smaller masses as a result of the present study. The 
comparisons made in this work will help us to calibrate 
different models of angular momentum evolution. 

5. First comparisons among models 

5.1. The role of the D-burning 

An uncertainty that deserves a separate discussion is the 
abundance originally assigned to deuterium. Deuterium 
burning, at least in the classic scenario, takes place in the 
early phases of the pre-MS evolution, when the central 
temperatures reach T~10^K. A higher deuterium abun- 
dance does not alter the track location, but stretches the 

^ Fig.Qand all others that use theoretical results were made 
with rotating models 



duration of the D-burning phase itself, thus determin- 
ing a delay in the evolution. We calculated for compar- 
ison a set of tracks with double deuterium abundance, 
A'(Z))~4x 10~^, that is the highest value acceptable on 
observational grounds (Linsky 1998). In comparing the 
two sets of tracks, we note the following: 

1. the track location does not change; 

2. for Af>0.5Afo, the high-D evolutions are slower, but 
the age difference is in all cases shorter than 200,000yr; 

3. for M<O.AMq, the age differences between tracks 
of the same mass are no longer negligible, from 
-250,000yr for M=O.AMq to -600,000yr for 
M^0.2Mq. 

In terms of the general conclusions found during the anal- 
ysis of the ONC stellar population, increasing the D abun- 
dance would have no effect on the mass distribution, while 
the age distribution would be shifted to slightly older ages. 

5.2. Role of non-grayness 

The gray tracks of models with a~1.0 are systematically 
hotter than their non-gray counterparts, although the dif- 
ference in Toff varies with the mass. For the lowest masses 
of our sample, i.e. M<0.2Mq, the high gravities allow the 
differences to remain within ATeff~100K. This difference 
increases to ^250 K at 0.6Mq and reaches a maximum 
of ^400 K in IMq models only at the end of the Hayashi 
track. The stars for which the pre-MS tracks are most sen- 
sitive to the boundary conditions are those with masses 
in the range O.AMq < M < IMq: the difference in T^e 
slightly increases along each track, and reaches a maxi- 
mum of ATcff ~ 400K when the radiative core is formed. 
The Tcff of more massive stars are less influenced by the 
atmospheric treatment. These differences in Tes lead to 
the assignment of a different mass to a given observed 
star. We verify that, consistent with the above discussion, 
this effect is negligible for the lowest masses, while it leads 
to differences of the order of -O.IMq for M~0.3-0.5Mq, 
where, as we shall see, the bulk of the ONC population is 
found. For larger masses the differences are larger because 
the tracks for different masses are closer to each other. 

5.3. The convection model 

The tracks corresponding to the three values of a used 
in the present investigation are shown in the left panel 
of the Fig. n The observed stellar loci are also reported. 
In the early phases, with the exception of the D-burning 
phase, gravitational contraction is the only source of en- 
ergy. The densities and temperatures increase until the 
central regions become radiatively stable; shortly after, H- 
burning takes over as the main energy source. We see from 
the left panel of the Fig. ^that the a2.0 and a2.2 tracks 
are systematically hotter, the differences being larger for 
higher masses and smaller for older ages. This behavior 
can be understood on the basis of the different degree of 
overadiabaticity present in the various masses at different 
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Fig. 2. Left: Variation with luminosity of the overadiabaticity of the stellar layer at t—10 of some pre-MS models 
calculated with different values of the parameter a entering the mixing length expression. Right: The variation with 
luminosity of the fraction of the mass of the star where (V — Vad) > 10~*, for the same masses shown in the left 
panel. For clarity we show only models al.O and a2.0. 



ages. The exact temperature profile depends on the overa- 
diabaticity (V — Vad), which is defined as the excess of 
the effective temperature gradient in comparison with the 
adiabatic temperature gradient, and it is only noticeably 
different from zero at the border of the convective zone, 
where convection becomes inefficient. The two panels of 
Fig. 121 show, respectively, the evolution of the overadia- 
baticity at T=10 and of the width (in solar masses) of the 
external region of the star where (V — Vad) > 10"''. For 
clarity we report only the al.O and a2.0 models. The lu- 
minosity is on the abscissa as a time indicator. A detailed 
inspection of Fig. |21 shows: 

1. the overadiabaticity at r=10 of the al.O models is sys- 
tematically higher than their a2.0 counterparts. This 
can be understood on the basis of the intrinsic lower 
efficiency of the convective model adopted, due to a 
lower mixing length; 

2. due to the higher internal densities of the less mas- 
sive models (hence, greater convective efficiencies) , the 
overadiabaticity differences increase with the mass; 

3. in the less massive models the differences above tend 
to narrow with age; 

4. the extension in mass of the overadiabatic region is also 
systematically higher in the al.O models, and tends to 
shrink with age. 

We verified that models of the same mass belonging to 
the two sets of tracks follow the same L{t) relation. The 
differences in the location of the tracks are, therefore, to 
be totally ascribed to differences in the effective temper- 
atures. Since the interior of these structures is practically 
adiabatic (in the center, (V — Vad)<10~'' in all cases), the 
radius is mainly determined by the degree of the overadia- 



baticity. This explains why larger differences are found for 
higher masses. In the lowest masses the efficiency of con- 
vection increases at older ages, so that the sensitivity to 
the adopted model for convection is strongly reduced. The 
low mass tracks approach each other at low luminosities 
(see Fig.[T}. 

5.4. The lithium depletion 

Table [21 shows the lithium concentrations for our three 
sets of tracks and their non-rotating counterparts. We re- 
port the values for 0.7, 0.8, 0.9, 1.0 and 1.2 M© at lO^yr. 
A comparison between log N [Li] found with rotating and 
non-rotating models, keeping a fixed, is qualitatively in 
agreement with the abundances found by Mendes et al. 
1)1999(1 . i.e., rotating models provide greater lithium deple- 
tion especially for low-mass stars at the age in question. 
Another comparison between the abundances found with 
different values of a, keeping the rotation status fixed, 

Table 2. Lithium abundances at the age of lO^yr for 0.7, 
0.8, 0.9, 1.0 and 1.2 Mq. We report the values obtained 
with different sets of our non-gray models, with ("rot") 
and without ( "non-rot" ) rotation. The initial lithium con- 
centration is logN[Li] = 3.31. 



Models 




O.7M0 


O.SMq 


0.9Mq 


I.OMq 


1.2M0 


rot Q=1.0 




0.750 


2.607 


3.087 


3.226 


3.291 


rot a=2.0 




-0.781 


1.669 


2.267 


2.720 


3.139 


rot a=2.2 




-0.961 


1.466 


2.109 


2.592 


3.094 


non-rot a= 


1.0 


0.950 


2.650 


3.097 


3.228 


3.291 


non-rot a= 


2.0 


-0.524 


1.715 


2.295 


2.733 


3.143 


non-rot a= 


2.2 


-0.521 


1.712 


2.297 


2.733 


3.143 
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Fig. 3. Mass (left) and age (right) histogram based, respectively, on the masses and ages assigned to the observed 
stars in the ONC using three different sets of tracks, calculated with three prescriptions for convection, a— 1.0, 2.0 
and 2.2. 



confirms the results by D'Antona & Montalban H2U(J3fl : the 
lithium depletion of the HCE models is too large to be con- 
sistent with the lithium depletion observed by Soderblom 
et al. H1993|l and Garcia Lopez et al. H1994|l in young open 
clusters, which can be reproduced only by the LCE al.O 
models. On the contrary, the solar radius is reproduced 
only by the a2.0 model, and 2D hydrodynamic computa- 
tions indicate HCE in the pre-MS. We regard this result as 
an indication that the efficiency of convection in the pre- 
MS might be affected by other parameter (s). The "second 
parameter" affecting lithium depletion is identified as the 
stellar rotation rate by Siess & Livio ( f 997) and in the pa- 
pers by Ventura et al. I|1998(l and D'Antona et al. fSOOO). 
The first authors propose that a is smaller in fast rotat- 
ing pre-MS stars due to the twisting of convective cells, 
the others show that the action of the dynamo-induced 
magnetic field due to the interaction of rotation and con- 
vection modifies the structure of the convective layers and 
reduces lithium depletion. 

6. Comparison with the ONC stars 

6.1. Data from the literature 

In order to study angular momentum evolution in pre- 
MS phase, we will compare our sets of evolutionary tracks 
with observational data of the ONC stars. To accomplish 
this goal we need some key parameters, such as effective 
temperatures and luminosities (to infer masses and ages of 
the stars) and also the rotation period and an index that 
allows us to distinguish between different kinds of angular 
momentum evolutions. 

The ONC data we used have been kindly provided by 
Dr. Keivan Stassun, who has widely worked on the rota- 



tional properties of ONC (Stassun et al. HHMI EHil) • Our 

final study sample is composed of a combination of data 
from the following sources: 

1. Rotation periods: Stassun et al. 1)1999(1 . Herbst et al. 
(20021 : 

2. Effective temperatures and luminosity: Hillenbrand 
CTMTjl : 

3. Infrared continuum excess, A[I— K]: Hillenbrand et al. 
ijTOMjl : 

4. Call equivalent width, EW(Call): Hillenbrand 
l(TM7jl . Hillenbrand et al. ifTU^ : 

5. X-ray luminosities: Stassun et al. I|2()()4|l . 

The rotation period diagnostic was obtained by photome- 
try, interpreting as rotation the periodic photometric vari- 
ability, probably due to the presence of stellar spots. 

The IR excess (A[I— K]) was obtained from combined 
optical and infrared photometric data. Extinction and 
spectral typing uncertainties are the main sources of er- 
rors in A[I— K] estimates. For earlier spectral type stars 
( ^ K2) the errors are negligible (<0.05 mag), for spec- 
tral types in the range K2— M3, typical errors are between 
0.1—0.3 mag, and it is largest for the latest spectral types, 
where mis-classification causes relatively larger errors. 

The EW(Call) was obtained from the optical spec- 
troscopic study of Hillenbrand H1997|l and analyzed by 
Hillenbrand et al. H1998|l . Their measurements uncertainty 
is estimated at O.SA based on measurements of multiple 
spectra of the same star. 

In order to obtain more reliable values for Lx 
of ONC stars, Stassun i2QMi reanalyzed all archival 
Chandra/ ACIS observations of these objects using up- 
dated calibrations and including time-filtering of flares. 
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Fig. 4. The comparison between the age distributions of the observed stars with assigned mass lower and higher than 
Mtr. The Mtr value is O.SMq for LCE models and O.SSMq for HCE models. 



The effective temperatures and luminosities were ob- 
tained from optical spectroscopy and photometry. 

6.2. Derivation of masses and ages 

For each of the three sets of tracks previously discussed, 
we assigned to each observed point a mass and an age by 
linearly interpolating between the two nearest tracks. An 
inspection of Fig. ^ can help understand, at least qual- 
itatively, the differences that we should find by varying 
a. For each observed star in the a2.0 and a2.2 sets, we 
assign a systematically smaller mass than in the al.O mod- 
els, hence a younger age (we recall that this evolutionary 
phase is governed by gravitational contraction, that pro- 
ceeds on a Kelvin-Helmoltz time scale tkhoc^). If the 
a2.0 or a2.2 sets are used, we therefore expect a mass dis- 
tribution shifted to lower masses and, on the average, a 
younger population. Typical internal errors are ^ 0.2 dex 
in log(L/L0) for all spectral types and are ^ 0.02 dex in 
logTcff for late-type (K— M) stars, but increase towards 
earlier spectral types. This leads to an uncertainty in the 
determination of mass that is ^ O.IMq for M<0.5Mq, 
and gradually increases to ~0.2Mq for M^IMq. The at- 
tribution of age is mainly influenced by the uncertainty on 
the luminosity, that makes the age uncertain by '^IMyr 
at the age of 1 Myr. This poses the problem of whether 
the age distribution we find should be considered either 
as the result of a burst of star formation or as a real in- 
dication of age differences from star to star. In the course 
of the investigation we favour a statistical interpretation 
of data as an indication of some age evolution, based on 
the evolution of rotation periods. 

The left panel of the Fig. O shows the mass distribu- 
tion of the observed stars, obtained by using the three sets. 
Only masses M<1.6 AIq are plotted because they repre- 
sent most of the stars in the sample. The mass function 
for the al.O set peaks in the mass interval 0.3— 0.4Mq, 
but we also note the presence of a significant group of 
stars with masses in the range 0.2— 0.3Mq and another 



group in the interval of 0.6— 0.9 Mq. For the HCE models 
this latter group of objects becomes less relevant and the 
mass function peaks in the mass interval 0.2— O.4M0. The 
right panel of Fig.|31confirms that the age distribution de- 
pends on the choice of a. We note that a slightly younger 
population is obtained as the value of a increases, and, 
in the q;2.0 and q;2.2 cases, a very young group of stars 
appears at ages ^1— 2xl0^yr, but is not present for the 
LCE (al.O) models. This can be understood by consid- 
ering, in the right panel of Fig. ^ the relative location of 
observed points and theoretical isochrones. The age differ- 
ences are also due to the different slope of the isochrones 
corresponding to the two sets of tracks, which, in turn, are 
related to the differential variation of temperature with 
mass, already outlined in tl5.3l In the HCE models there 
is a group of objects with ages clustering around 100, 000 
yr . For the a2.0 case this young group is made up of ~40 
stars and for the a2.2 set its presence is still more evident, 
exhibiting ~70 objects. 



In ^7.1\ we will define a transition mass Mtr for HCE 
and LCE models, on the basis of the rotational periods 
distribution. Here we compare, in Fig. 01 the age distribu- 
tion for two different ranges of mass, M>Mtr and M<Mtr- 
We see from the left panel of Fig. 21 that in the al.O set 
the two populations show a similar distribution. On the 
contrary, the age distribution of the two groups for HCE 
models is very different. Thus, while the existence of a 
group of younger stars would be possible in the formation 
history of the ONC, it should be present for any mass. The 
discrepancy in the age distribution may be again an indi- 
cation that the al.O set provides a better description of 
the ONC stellar population, in agreement with the lithium 
depletion discussed in ii5.4l and with the previous analysis 
made by D'Antona & Montalban (j2003|l . In any case, this 
certainly is not final and we proceed with the analysis by 
using the three sets of tracks. 
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Fig. 5. The period histogram of all observed ONC objects. 
The total distribution of periods has a bimodal character, 
showing a primary peak of fast rotators (P'--^2 days) and 
a secondary peak of slow rotators (P~8 days). 



6.3. Comparison with gray models 

We compare now the effects introduced by the two major 
factors, i.e. convection and boundary conditions. Although 
for masses M>1Mq the treatment of convection is more 
relevant than the boundary conditions adopted in deter- 
mining the effective temperature of the tracks (Montalban 
et al. l2()()4|l . for the interesting range of mass for the ONC 
comparison, namely 0.2— 0.4 Af©, and the non-gray mod- 
els, convection is important mainly in the early evolution- 
ary phases. At later phases (ages >lMyr), convection be- 
comes more adiabatic and the non-grayness becomes the 
main factor affecting the track location. 

For gray atmospheric treatment, keeping a fixed to 
1.0 would concentrate the same mass distribution at 
M~0.2— 0.3 Mq. There would also be a considerable re- 
duction of the population with masses in the range 
Q.6<M/Mq<1. a larger effect would be obtained if we 
had used gray models with q;=2.0. In this case, the mass 
function would be peaked in the range 0.1<M/Mq<0.3, 
and the average age of the observed stars would be slightly 
younger than ~lMyr. This latter result was obtained by 
Herbst et al. H2002fl during their analysis of the ONC pop- 
ulation using the tracks by D'Antona & Mazzitelli (1994) 
that use gray approximation and the very efficient FST 
model for convection. 



Table 3. Main physical parameters of the present gray 
(G) and non-gray (NG) models. A/'< (A/'>) is the percent- 
age of the Nt stars that have mass less (greater) than Mji- 
for different rotation periods (P). See text for details. 



Models 


Mass 
range 

(Mo) 


Age 
range 

(Myr) 


Mtr 

(Mo) 


P<4d 


P>6d 


P>6d 


G Q=1.0 


0.2-0.3 


0.6-2.5 


0.35 


65% 


19% 


53% 


G Q=1.5 


0.1-0.3 


0.3-1.3 


0.25 


65% 


19% 


54% 


NG a=1.0 


0.2-0.4 


1-2 


0.5 


63% 


18% 


53% 


NG a=2.0 


0.2-0.4 


0.6-2.5 


0.35 


67% 


18% 


55% 


NG a=2.2 


0.2-0.4 


0.4-1.6 


0.35 


67% 


23% 


54% 



7. Stellar rotation in the ONC 

7.1. The dichotomy in period distribution for different 
mass ranges 

We examine in Fig. the distribution of the observed ro- 
tation periods. We note the presence of a primary peak 
corresponding to fast rotators with 1 ^ P(d) ^ 3, and a 
secondary peak at P~8d. The former can be associated 
with spin up due to the conservation of the total angular 
momentum. The latter indicates the presence of a mech- 
anism acting to prevent stellar spinning up, at least in 
the early evolutionary phases. Attridge & Herbst H1992|l . 
Choi & Herbst Ijl996|) suggested that this can be due to a 
"disk locking" , due to magnetic coupling between the star 
and the disk (K6nigl[TMT|). Following Herbst et al. ll^UH^ 
we investigate in detail the rotational status of the vari- 
ous masses involved. The rotational properties of the stars 
vary considerably with mass: stars with masses larger than 
a value Mtr have a clearly bimodal distribution, while the 
less massive sample contains only a tail of slow rotators. 
This behavior was first observed by Attridge & Herbst 
l(TM^ and discussed by Herbst et al. ^UU^. We define the 
"transition mass" (Mti), which depends on the track set 
chosen for the analysis, in order to maximize the effect of 
the bimodality. Although the dichotomy does not depend 
on the chosen set, the transition mass varies according to 
it. For LCE models it is O.SM©, while a reasonable value 
is O.35M0 for HCE models. If we had used our gray mod- 
els with a=1.5, Mtr would have been still smaller, namely, 
O.25M0, in agreement with the findings by Herbst et al. 
(12002), who used the HCE FST models by D'Antona & 
Mazzitelli p994|l . In Fig. we show the histogram of pe- 
riods, respectively, for stars less and more massive than 
Mtr- The secondary peak at P'^ 8d, already seen in the 
Fig. 13 is present only in the population at M>Mtr, while 
the low mass objects show a clear trend towards short pe- 
riods. Table 01 shows the percentages of slow and fast ro- 
tators (here defined by the limitations P>6d and P<4d). 
Fast rotators contains more than 60% of masses <Mtr. On 
average, only 20% of the masses <Mtr and ^54% of the 
masses >Mtr have P>6d. 

This dichotomy indicates that either (i) disk locking is 
responsible for the presence of the secondary peak, and 
stars with M>Mtr tend to be embedded in their disk 
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Fig. 6. Period iiistograms showing the dependence on mass of the period distribution of the ONC objects. Stars more 
massive than Mtr have a bimodal period distribution and their less massive counterparts rotate faster and exhibit a 
unimodal distribution. This behavior can be seen independently of the choice of a. The Mtr value is the same used in 
Fig.H 
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Fig. 7. The mass distribution of the sources which according to our analysis rotate fast since early evolutionary phases 
(solid lines) compared to those that are slow rotators in the current epoch (dotted lines). This comparison is shown 
for the three sets of models (al.O, a2.0 and a2.2). 



longer than their low mass counterparts, (ii) the locking 
time is similar, but the masses >Mtr evolve faster and a 
larger fraction of their pre-MS lifetime is locked, or (iii) 
the "locking period" of the group with M<Mtr is signifi- 
cantly lower than ~8d. Our analysis confirms possible in- 
terpretations of the observed distribution of periods given 
by Herbst et al. j2002) (see their Fig. 15). The uncertainty 
on the convection model simply alters the transition mass 
from ~O.25M0 to a maximum of 0.5M© for LCE models. 

7.2. Disk locking and the disk lifetime 

Following the suggestion by Herbst et al. I|2U(J2II , that the 
longer period peak in the distribution indicates that some 
stars are locked in their disks with a period near 8 days, 
we simply considered the stars with periods larger than a 
threshold period (Pthrosh)j which we put at 8d, as locked. 
For stars at P<Pthresh, unlocked stars according to our cri- 
terion, we determine the epoch at which their period was 



equal to 8 days. This would be the time at which the stars 
would have lost their disks, and began a constant angular 
momentum evolution. Here we consider that angular mo- 
mentum losses by magnetic braking are negligible at the 
pre-MS, since its timescale is much longer than the evolu- 
tionary timescale during pre~MS evolution. The temporal 
variations of radius and angular velocity were determined 
on the basis of our tracks once the mass was assigned. 
Following this hypothesis, we found some stars which had 
P=Pthicsh at an age younger than lO^yr. In our interpre- 
tation, these stars have lost their disks very early and can 
be considered to have evolved without a disk. In this way, 
we identify three distinct populations: 

1. early fast rotators - stars locked only for ages < lO^yr, 

2. slow rotators - stars probably still disk embedded, 

3. moderate rotators. 

Stars in the last group may have lost their disks at ages 
greater than lO^yr. They represent, on average, ^--^45% of 
the stars of the whole sample and this percentage do not 
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early fast rotators 
slow rotators 





0.5 1 1.5 0.5 1 1.5 0.5 1 1.5 

M/Mj, M/Mj, M/Mg 

Fig. 8. The observed infrared excess of our sample stars plotted against their inferred mass according to the three 
sets of models used in this work. 



change significantly with the choice of a. As long as the 
assumption that the locking period is independent of mass 
is valid, the disks seem to survive longer for higher masses. 
For all sets of tracks, the percentage of early fast rotators 
is --40% for M<O.4M0, and gradually drops below -12% 
for M>O.6M0. 

The mass distribution of the slow and early fast rota- 
tors is shown in Fig. [T] We note, in particular, the max- 
imum at 0.6-0.8 Mq (for LCE models) and 0.3-0.5 Mq 
(for HCE models) characterizing the slow rotator popula- 
tion (dotted line). We can compare the fraction of early 
fast rotators identified in the ONC with the non-accreting 
fraction of TT and Brown Dwarfs in star forming regions 
of similar ages (p Oph and Taurus, Mohanty et al. I2005'il . 
The lower limit to this non-accreting fraction is ^---35%, 
not very different from our result. 

In our analysis we use the rotation period as an indi- 
cator of the presence of a disk surrounding these stars. In 
order to test the reliability of this hypothesis, we should 
use several observational indicators of the presence of disk 
and accretion, like the infrared excess A[I— K], the equiv- 
alent width of Call line, the excess in the L-band and 
Hq,, Ca, O emission lines. Near IR excess as disk indica- 
tors and EW Call as accretion diagnostic must be used 
with caution f Hillenbrand II 997|l . But, as this is mainly a 
theoretical work, we will check what is already in the lit- 
erature for ONC stars, namely IR excess A[I— K] and EW 
Call, mainly as additional arguments. 

It is expected that still locked stars have A[I— K]>0.3, 
and those that evolved without disk should have infrared 
excess significantly lower than this threshold value (Herbst 
et al. l2()()2|l . We report in Fig.|Hlthe observed stars on the 
plane A[I— K] vs. mass. We can see that sources that we 
identified as still locked (slow rotators - crosses) are mainly 
concentrated above the A[I— K]=0.3 line, while those that 
evolved without a disk (early fast rotators - open squares) 
lie mainly below it, for the three sets of models. This 
straight correlation between the infrared excess and our 
derivation of the presence of a disk agrees with our theo- 
retical considerations. 



The equivalent width of Call lines is commonly used 
as an indicator of an active accretion process. For accret- 
ing stars we expect emission lines and EW(CaII)< — 1, 
while for non-accreting objects we have EW(CaII)>l 
(Flaccomio et al. I2003b|l . We cannot expect a direct cor- 
relation between A[I— K] and EW(Call), because some 
stars might still have a disk, although no longer ac- 
creting. Yet, we expect that the observed stars with 
EW(CaII)<-l and A[I-K]>0.3 should have a disk sur- 
rounding them, and the disk locking mechanism should be 
active. We could identify ~40 stars with known rotational 
periods that satisfy both requests, about 30% of which are 
actually identified as still disk embedded according to our 
criterion (P>Pthrcsh=8d). Part of the remaining have long 
rotational periods, very close to Pthicsh, which suggest the 
presence of a disk. 

7.3. An alternative view: the role of the magnetic field 

The idea that disk-locking is responsible for the dichotomy 
observed, which eventually leads to the presence of two 
peaks in the distribution of periods of the ONC stars, was 
recently criticized by Barnes H2003|l , who argued that the 
possible role of disks can be only to set an initial distri- 
bution of periods, since disk-locking should affect all the 
stars equally. 

Based on the observed color-period diagrams of sev- 
eral open clusters (e.g. IC 2391, a Per, Pleiades, M34, 
etc.), Barnes ( 2003) found that a double population with 
distinct rotational properties characterizes any stellar as- 
sociation. Further, he found that a systematic trend with 
age is apparent, namely that older clusters have a smaller 
number of rapid rotators, that eventually disappear at 
ages ~800 Myr. 

Barnes H2003|l interpreted these observations as a re- 
sult of a different morphology of the magnetic field config- 
uration. According to his suggestion, rapid rotators have 
small scale magnetic fields associated with their convec- 
tive region, which cannot be anchored either to the in- 
ner radiative core or to the star's external layers. All 
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Fig. 9. The temporal evolution of the periods of our three sets of models with mass in the range 0.2 Mq<M<OA Mq 
evolved starting with an initial angular momentum calculated on the basis of the Kawaler (jl987|l prescription (solid 
lines), and with the same values multiplied by a factor of 3 (dashed lines). The stars that we suppose to be evolved 
with a constant angular momentum, i.e., without a disk, are also shown. The o symbols identify the mass range 
O.2<M/M0<O.3, ■ is used for 0.3<Af/MQ<0.4, and x for all the other masses. 



the fully convective stars should belong to this group. 
Conversely, the slowly rotating stars are characterized by 
large scale magnetic fields, probably associated with the 
presence of an interface dynamo between the external con- 
vective region and the internal radiative zone. In this case 
the process of spinning down the star is much more effi- 
cient. The dynamo, that is probably created by the decou- 
pling between the convective and radiative zones, anchors 
the spun-down convective envelope to the rapidly rotat- 
ing core, thus favouring a constant migration of the stars 
belonging to the rapidly rotating group to the slowly ro- 
tating sample. This should explain the complete absence 
of the fast rotator sequence at old ages (see Barnes BDU51 
Fig. 1). Due to the young age of the ONC, both popula- 
tions should be present there, as we find. This interpreta- 
tion is not correct if we take our results at face value: if we 
rely on our attribution of masses and ages, the great ma- 
jority of the stars observed are indeed fully convective, i.e. 
should all belong to the rapidly rotating sequence. We sug- 
gest that the magnetic field itself plays a role in inhibiting 
convection (Cough & Tavler lTM^ Moss ll9b'5l Ventura et 
al. 11998(1 and favours an earlier appearance of a radiative 
core in the stars having M>Mtr. Consequently, we do not 
reject the idea that the presence of a double population 
is indeed due to the magnetic field configuration, rather 
than to the effects of a disk-locking mechanism. 

7.4. A constant angular momentum evolution? 

For those stars that we identified to have lost the disk in 
their early evolutionary phases, i.e., the early fast rotators, 
we may test the hypothesis that they evolved at constant 
angular momentum from the beginning. The distribution 
of these objects in the P vs. inferred age plane is shown 
in Fig. El for a=1.0, a=2.0 and a=2.2 models. We used 
different symbols as mass identifiers. We can see a clear 



trend towards shorter periods for older ages, especially for 
the al.O models, possibly indicating angular momentum 
conservation. We checked the possibility of reproducing 
the observed rotational pattern with age by means of our 
rotating models. This approach allows us to find out the 
range of initial angular momenta that, for each mass, must 
be used to calculate the models. We limit this analysis to 
the subsample of the mass range given in Table 13 We 
divided the observed sources into three classes of mass, 
indicated in Fig. with open circles masses in the inter- 
val O.2<M/M0<O.3, with full squares O.3<M/M0<O.4 
and with crosses all the remaining. The solid lines indi- 
cate the temporal variation of the rotational periods ac- 
cording to the evolution of our three sets of models with 
masses M=0.2, 0.3, 0.4 Mq calculated by assuming an ini- 
tial angular momentum following the prescriptions given 
in Eq. We note that the temporal evolution of the pe- 
riods vary with the parameter a, as it affects the radius 
of the stars. These curves can only reproduce the upper 
envelope of the observed loci, but, particularly at the ages 
shared by the bulk of the ONC stellar population, they 
lead to rotational velocities too slow with respect to most 
of the observed values. 

To reproduce the rotation period of the fastest stars, 
we need to use an initial angular momentum at least three 
times larger than that prescribed by Kawaler l|1987|l . if we 
use LCE models, and even larger for HCE models. To fully 
bracket the observed periods it is necessary to assume a 
distribution of initial angular momenta Jjn, at least, in the 
range Jkaw<>/in<3 Jkaw This result can be used to extend 
the Kawaler H1987|l prescription to the very low mass stars. 

7.5. The X-ray emission of the ONC stars 

Flaccomio et al. H2003al I2003b|l and Stassun et al. pUIMjl 
performed deep analyses of the archival Chandra data and 



Landin et al.: Non-gray rotating stellar models and the evolution in the ONC 



15 




Fig. 10. Top: Lx luminosity (ergs/s) plotted against mass (obtained with our three sets of models) for the ONC 
observed sources with known rotational periods. Bottom: The same as top, with bolometric luminosity on the vertical 
axis. 



derived the Lx luminosity of all the sources included in 
the Hillenbrand |1997J sample. Their goal was to correlate 
Lx of pre-MS stars with the factors most likely driving the 
X-ray emission itself, i.e., accretion and rotation. Their 
main finding was the lack of a clear correlation between 
Lx and rotational period. They interpreted this result as 
evidence that the ONC pre-MS stars are indeed in the 
"super-saturated" regime of the rotation-activity relation- 
ship. This seems to be confirmed by the average value of 
the fractional X-ray luminosity log(Lx/iboi)'^— 3.6, that 
is slightly smaller than the main sequence saturation value 
of log(Lx/Lboi)^-3. 

Concerning the relationship between accretion and X- 
ray emission, Stassun et al. (2004) found that accreting 
stars have X-ray luminosities on average lower than their 
non-accreting counterparts, as a possible result of X-ray 
extinction by circumstellar gas in magnetospheric accre- 
tion columns. We could not find any clear correlation be- 
tween EW(Call) and Lx, although this might be a con- 
sequence of the smaller sample of stars with evidence of 
accretion (i.e. EW(CaII)< — 1) used in the present work, 
compared to the complete sample by Hillenbrand H1997|l 
and Hillenbrand et al. H1998|l analyzed by Stassun et al. 

pro . 

We used our determinations of mass to look for any 
relationship between stellar mass and X-ray luminosity. 
In agreement with Flaccomio et al. H2003a|l . we find that 
Lx is correlated with mass (upper panels of Fig. ITIUl . 
This trend is not due to a qualitative difference in the 
X-ray emissions, but rather to a general correlation be- 
tween mass and bolometric luminosity (lower panels of 
Fig. I10|l . This is confirmed in Fig. ^] where we see that 
the Lx/Lboi ratio is practically independent of mass for 
the three a models we have used. We note the high dis- 
persion around the average value of log(Lx/Lboi)'^— 3.6, 
present at the lowest masses, that is probably connected 
to the lower luminosities of these objects. 



8. Conclusions 

We use the observed stellar population of the ONC to 
test our pre-MS evolutionary tracks, to better understand 
the main physical properties characterizing the evolution 
of young stars. By comparing the location of the tracks in 
the HR diagram with the position of the observed objects, 
we assign to any single star a mass and an age, for three 
different convection efficient models. 

We find the well known result that the treatment of 
convection is generally the most relevant physical input 
in determining the T^s of pre-MS tracks in the HR dia- 
gram. The boundary conditions adopted also play a non- 
negligible role in determining the path followed by the evo- 
lutionary sequences on the HR plane. Gray models are sys- 
tematically hotter than their non-gray counterparts. This 
efi'ect, for the range of masses and ages at which most 
of the ONC population is found, has a similar quantita- 
tive effect as a change of the convective model. The use of 
non-gray models is recommended to describe these early 
evolutionary phases. 

On the observational side, we find that the bulk of 
the observed stars in the ONC have masses in the range 
0.2Mq<M<0.4Mq, for all non-gray models. The age dis- 
tributions are more affected by the choice of the MLT 
parameter a. Ages are 1—2 Myr from the al.O set, 
0.6-2. 5Myr for a2.0 set and 0.4-1.6Myr for q;2.2 models. 
This study confirms the presence of a dichotomy in the 
rotational properties between the objects with M<Mtr, 
whose period distribution peaks at short values, and stars 
with M>Mtr, that present a secondary peak at P~8d. The 
transition value of mass between the two populations is at 
M~O.5M0 for LCE and at M~0.35Mo for HCE models. If 
disk-locking is responsible for the secondary peak observed 
in the overall period distribution, these results can be in- 
terpreted by assuming either that the masses M<Mtr lose 
their disk earlier, or that their locking period is shorter. 
The X-ray emission shows no correlation with period, sup- 
porting the suggestion that these stars are indeed in the 
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Fig. 11. The fractional X-ray luminosity as a function of mass for the same stars and models shown in Fig. 1101 



super-saturated regime of the rotation-X ray luminosity 
relationship. The correlation of the X-ray flux with mass 
appears to be the consequence of the increased average 
luminosity of more massive objects. 

For the low-mass stars that presumably evolved with- 
out a disk, we find that our results are consistent with 
an evolution-conserving angular momentum. The compar- 
ison between the model period evolution and the observed 
values suggests that initial angular momenta at least 3 
times larger than those found by means of the Kawaler 
H1987|l law are needed for stars with mass in the range 
O.2M0<Af<O.4 A/q. This analysis was not possible for 
the higher masses sample, due to the lack of a statistically 
meaningful sample of rapidly rotating stars. 

The idea that the double population of the ONC can 
be explained on the basis of a different morphology of the 
stellar magnetic fields seems to be ruled out by the fact 
that almost all the observed sources are fully convective, 
according to our interpretation. However, a mechanism 
(e.g. the magnetic field itself) inhibiting convection might 
favour an earlier appearance of the radiative core, at least 
in some of the stars. In our analysis we found other in- 
dications that convection in the pre-MS may be affected 
by other parameter(s): although 2D hydrodynamic simu- 
lations predict HCE in the pre-MS, we find two results in 
favour of LCE: (1) the lithium depletion in HCE models 
is too large to be consistent with the pre-MS depletion 
shown in young open clusters; (2) the age distribution de- 
rived from HCE models for two groups of smaller and 
higher masses is very different. It may be that the lower 
convection efficiency needed in the pre-MS is due to the 
structural effect of the dynamo-induced magnetic field, as 
suggested by Ventura et al. H1998|l and D'Antona et al. 
H2()()()|l . This possibility leads us not to dismiss the idea 
of an earlier appearance of a radiative core in the M>Mtr 
group. 

The evolutionary tracks (from 0.085 to I.6M0) and 
isochrones (from 2x10^ to IxlO^Myr) are available from 
the following web site: 

www.mporzio.astro.it /^tsa 
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